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ABSTRACT 

We discuss the linear response to low-frequency tidal forcing of fluid bodies that are 
slowly and uniformly rotating, are neutrally stratified and may contain a solid or fluid 
core. This problem may be regarded as a simplified model of astrophysical tides in 
convective regions of stars and giant planets. The response can be separated into non- 
wavelike and wavelike parts, where the former is related instantaneously to the tidal 
potential and the latter may involve resonances or other singularities. The imaginary 
part of the potential Love number of the body, which is directly related to the rates 
of energy and angular momentum exchange in the tidal interaction and to the rate 
of dissipation of energy, may have a complicated dependence on the tidal frequency. 
However, a certain frequency-average of this quantity is independent of the dissipative 
properties of the fluid and can be determined by means of an impulse calculation. 
The result is a strongly increasing function of the size of the core when the tidal 
potential is a sectoral harmonic, especially when the body is not strongly centrally 
condensed. However, the same is not true for tesseral harmonics, which receive a 
richer response and may therefore be important in determining tidal evolution even 
though they are usually subdominant in the expansion of the tidal potential. We also 
discuss analytically the low-frequency response of a slowly rotating homogeneous fluid 
body to tidal potentials proportional to spherical harmonics of degrees less than five. 
Tesseral harmonics of degrees greater than two, such as are present in the case of a 
spin-orbit misalignment, can resonate with inertial modes of the full sphere, leading 
to an enhanced tidal interaction. 

Key words: hydrodynamics - waves - planets and satellites: general - planets and 
satellites: interiors - planet-star interactions - binaries: close 



1 INTRODUCTION 

The tidal interaction of two bodies orbiting about their cen- 
tre of mass is one of the fundamental problems of theoretical 
astrophysics. It is relevant to a wide variety of systems, in- 
cluding close binary stars, the satellites of solar-system bod- 
ies, and stars orbiting the black holes in galactic centres. 
Interest in this subject has grown through the discovery in 
recent years of numerous planets in orbits of a few days or 
less around other starsQ It is likely that tidal interactions 
have affected the orbital and spin properties of many of these 
planets as well as providing a source of heat. Equally impor- 
tantly, they may also have led to the destruction of many 
planets that are not observed. 

Tidal theory involves a combination of celestial me- 
chanics and continuum mechanics. The celestial mechan- 



See exoplanet . eu or exoplanets.org. 



ics is non-trivial, and has been widely investigated us- 
ing highly simplistic descriptions or parametriz a tions of 



the fluid or solid behaviour fe.g. iDarwinl ll88(J: Ijeffrevs 


Il96l| 


IGoldreichl 19631: Goldreich & Soterl 1966: Alexander 


1973 


|]VIignardlll980l: lHutlll98ll: lEggleton, Kiseleva & Hut 


Il99a 


iMardling fc Linll2002l;lBarker fc Ogilviell2009l). In fact, 



the fluid dynamics of tidally forced bodies is a much 
richer and more difficult problem, even in a linear the- 
ory where the tide is treated as an inf initesimal d istur- 
bance. In the case of stars, the work of IZahnl l| 19771 . and 
references therein) identified two potentially important pro- 
cesses: the interaction of tidal bulges with turbulent con- 
vection and the excitation and da mping of internal grav- 
ity w a ves in radiative zones (s e e also ISavoniie fc P apaloizou 
I l983l: IGoldreich fc Nicholsonl Il989l: iTerquem et al.l Il998l ; 
boodman fc Dickson! Il998l : iBarker fc Ogilviel |2010J ). Rota- 
tion was neglected in this work in order to simplify the calcu- 
lations. Its effects have been considered bv lWitte fc Savoniiel 
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(1999) and lSavoniie fc Wittd (|2002h within the 'traditional 
approximation', which is usually applicable to the radiative 
zones of slowly rotating stars. 

More recent work has examined the effects of the full 



Coriolis force on tides in convective zones ( Ogilvie fc Li: 



Wul l2005alfB iPapaloizou fc Ivanovl 
I Ogilvie fc Linl 
lOgilviei 120091; 

also 



20051; IQgilvi* 
120071; iGoodman fc Lackne 
Rieutord fc Valdettarol 



2010 



1995 



Savoniie. Papaloizou fc^Ubertsl 

ISavoniie fc Papaloizoul 1 19971 : Papaloiz ou fc Savoniid 119971 ) 
All of these papers show that inertial waves, for which the 
Coriolis force provides the restoring force, can be excited 
by tidal forcing and may provide an important route for 
tidal dissipation in rotating bodies. However, even the 
linear theory of inertial waves is a very intricate problem 
and these authors come to different conclusions. In some 
cases the response is dominated by wave attractors, critical 
latitudes or other singularities, while in others global 
normal modes are excited. Different predictions are made 
for the frequency-de pendence of th e tidal dissipation rate . 
For example, while lOgilvie fc Linl i|2004l ). lOgilviei (|2009l ) 



and iRieutord fc Valdettarol (|201C ) allow inertial waves 
to undergo multiple reflections within a spherical shell 
before being damped by viscosity, and obtain a highly 
freque ncy-dependent dissipation rate, IGoodman fc Lackner] 
(2009) do not allow any reflections of waves after their 
generation from the boundary of the planetary (or stellar) 
core, but instead consider nonlinear damping and obtain 
a smooth dissipation rate that inc reases with th e fifth 
power of the core size. Nevertheless, lOgilviei (2009) found 
that the frequency-averaged dissipation rate with multiple 
reflections shows the same dependence on core size. 

One purpose of this paper is to show that, in spite of 
this complexity, some broader aspects of the tidal response of 
rotating bodies are robust and insensitive to the details that 
distinguish between these calculations. Even if the linear 
theory of inertial waves in a perfect spherical shell cannot 
be applied reliably to the turbulent convective zone of a star 
or planet, the results presented here on frequency-integrated 
responses should still be valid. 

A second purpose of this paper is to collect some ana- 
lytical results on the tidal response of a homogeneous fluid 
body. When potential components other than the usual 
I — m — 2 spherical harmonic are considered, a richer re- 
sponse is possible and numerous resonances can occur with 
large-scale inertial modes, especially in systems with a spin- 
orbit misalignment. Similar behaviour can be expected in 
more realistic models, and the astrophysical consequences 
require further investigation. 



v s — e r R — e r x (e r x VS), 
v t = -er x VT. 



(2) 
(3) 



This nomenclature is used, for example, in seismology 
(e.g. lLapwood fc Usamilll98ll ) and stellar oscillations (e.g. 
ISmeveral2010r i. and is related to the theory o f vector spher- 
ical harmonics (e.g. Morse & Fcshbach 1953). We have 



e r ■ v — R, 

e r ■ V X (— e r x v) 



-vis, 



e r ■ V x v 

where 



-VST, 



d0 



( s[n9 Yo 



+ 



r 2 sin < 



(4) 
(5) 
(6) 

(7) 



is the horizontal part of the Laplacian operator. On each 
sphere r — constant > 0, the operator VJ^ has a complete set 
of eigenfunctions (spherical harmonics) and can be inverted 
uniquely except for the addition of an arbitrary constant. 
Therefore, given a differentiable vector field v, equations 
([4]|-((6]) determine R, S and T uniquely except that arbitrary 
functio ns of r only ca n be a dded to S and T. 

In lOgilvie fc Linl (|2004 ) and lOgilviei (|2009h the veloc- 
ity field is separated into spheroidal and toroidal parts and 
expanded in spherical harmonics. The coefficients a n and 
b n refer to the spheroidal part (R and S), while th e coeffi - 
cients c„ refer to the toroidal part (T). In lOgilviei ([2009), 
for example, we have 



S = r 2 J2bu(r)Y^(0, 



where the sums are over integers n > m > 0, 



(8) 



(9) 



(10) 



(11) 



(12) 



and the real part of these expressions is to be taken after 
they are multiplied by e -1 "*. 



(2n+ l)(n-m)! 



1/2 



P™ (cos 6) e 1 



47r(n + to)! 

is a spherical harmonic with the normalization 
\Y™(8,(j>)\ 2 sinfldfld^ = 1, 



2 SPHEROIDAL AND TOROIDAL VECTOR 
FIELDS 

We begin with a definition that is important for the analysis 
in this paper. Let (r, 8, <f>) be spherical polar coordinates, 
and e r — r/r the radial unit vector. A general differentiable 
vector field v can be represented in the form 

v = e r R - e r x (e r x VS) - e r x VT, (1) 

where R, S and T are scalar functions of (r, 6, <f>). Thus v = 
v s + v t , where the spheroidal and toroidal parts of v are 



3 TIDAL RESPONSE OF A HOMOGENEOUS 
ROTATING BODY 

Perhaps the simplest tidal problem is the linear response of 
a homogeneous fluid sphere to a tidal potential that is har- 
mon ic in time and p roportional to a solid spherical harmonic 
(see I Thomson! Il863l . although his analysis is mainly for an 
elastic solid rather than a fluid). To include the effects of 
rotation, we sh ould properly consid er an ellipsoidal figure 
of equilibrium (|Chandrasekharlll969TK ti des in a Maclaurin 
spheroid were analysed by iBrvanl (|l889T ) , although the full 
ramifications of his work have not yet been explored. In the 
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limit of slow rotation, however, it is adequate to neglect cen- 
trifugal effects and consider a spherical fluid body subject 
to the Coriolis force. 

The tidal gravitational potential experienced by a body 
can be written as an interior multipole expansion in solid 
spherical harmonics of degrees I > 2; this series converges 
more rapidly when the tide-raising body is more distant. For 
each value of I, the response of a rotating body also depends 
on the order m of the spherical harmonic, which satisfies 
—l<m<l. (The response for negative values of m can be 
deduced from that for positive m by complex conjugation if 
the frequency is also changed in sign.) 

We collect in Appendix[5]some analytical results on the 
linear response of a homogeneous sphere to forcing by solid 
spherical harmonics with degrees I = {2,3,4}, i.e. interior 
quadrupolar, octupolar and hexadecapolar potentials. (No 
tidal force is generated by I — or Z = 1.) This analysis 
makes use of the decomposition of the tidal response into 
non-wavelike and wavelike parts, developed in Section 14.41 
below. 

In general, the tidal response involves resonances with 
inertial modes of the sphere when the tidal frequency 
matches the eigenfrequency of an appropriate inertial mode. 
The eigenfrequencies of inertial modes in the frame that ro- 
tates with the angular velocity Q of the body are dense in 
the interval (-2Q, 2fi) i|Greenspanlll968l) . However, the se- 
lection rules determined by the spatial structure of the iner- 
tial modes and the forcing functions imply that any spher- 
ical harmonic potential resonates with only a finite number 
of inertial modes. Indeed, as found in Appendix [J] the sec- 
toral harmonics (\m\ — I) do not resonate with any inertial 
modes, while the tesseral harmonics (\m\ < I) do have such 
resonances. 

However, the resonances of the quadrupolar tesseral 
harmonics Y 2 ° and are only formal resonances and should 
not lead to enhanced dissipation in linear theory. The for- 
mally resonant modes have zero frequency in an inertial 
frame of reference; they are present because of the existence 
of a continuous family of equilibrium solutions for uniformly 
rotating bodies that have angular velocity vectors that dif- 
fer in magnitude and/or direction. In the case m — the 
zero-frequency mode is the 'spin-up' mode, which consists 
simply of an infinitesimal uniform change in the angular ve- 
locity of the fluid. In the case m = 1 it is the 'spin-over' 
mode of frequency u> — — fl (i.e. zero frequency in the in- 
ertial frame), which consists of an infinitesimal tilt of the 
rotation axis of the fluid. Each mode involves a change in 
the exactly conserved angular momentum of the body and 
therefore has zero frequency. The formal resonances of Y° 
and Y% with these modes could not produce any tidal dis- 
sipation because the modes concerned involve uniform rota- 
tion and the resonance occurs at zero forcing frequency (in 
the inertial frame). It may be surprising that the spin-up 
mode, which involves a change in the angular velocity of the 
fluid, can be excited by an axisymmetric tidal potential pro- 
portional to y 2 °i which cannot change the axial component 
of angular momentum of the body. In fact, even in the ab- 
sence of a torque, the spin-up mode is excited because the 
axisymmetric tidal deformation modulates the moment of 
inertia of the body. 

The most important harmonic present in the tidal po- 
tential is usually considered to be Y? , and the fact that 



it cannot resonantly excite inertial normal modes in a 
homogeneous full sphere h as been noted previously (e.g. 
iGoodman fc Lacknerl hootf ). However, given that tesseral 
harmonics have a richer response in a homogeneous body, 
and therefore presumably also in an inhomogeneous body, 
other harmonics are worthy of consideration, even if they 
are subdominant in the expansion of the tidal potential. 

To illustrate this, we consider here the case of two bod- 
ies in a circular orbit of angular velocity (or mean motion) 
n. We consider the tide raised by body 2 in body 1, which 
rotates uniformly, and assume that the orbit is inclined with 
respect to th e equatorial plane of bod y 1. [ Rela t ed an alyses 
are given by iBarker fc Ogilviel i|2009l ) and [La] (|2012h .] We 
consider four spherical polar coordinate systems, all centred 
on body 1: 



Si 

S-2 
S 3 

S 4 



aligned with the orbit, rotating with the orbit 
aligned with the orbit, non-rotating 
aligned with the spin of body 1, non-rotating 
aligned with the spin of body 1, rotating with body 1 



In 5*1 the tide-raising body (as represented by its centre of 
mass) is stationary and lies in the plane 9 = n/2. The tidal 
potential is independent of time and can be expanded in 
solid spherical harmonics. Let us focus on terms of a fixed 
degree I > 2. Then harmonics of orders m = {I, I — 2, I — 
4, . . . , —1} are present, because symmetry about the plane 
6 = n/2 requires that I — m be even. In S2 the tide-raising 
body executes a circular orbit in the plane 8 = n/2. The 
tidal potential involves the same spherical harmonics but 
their angular frequencies are now mn. The rotation from 
S2 to S3 mixes all permissible values of m for a given I 
(as described by the Wigner D-matrix). In S3, therefore, all 
values of m with —l<m<l are present and the frequencies 
are {I, I — 2, I — 4, ... , —l}n for each m. Finally, in S4, all 
values of m are present and the frequencies are {1,1 — 2,1 — 
4, . . . , —l}n — mil, where fi is the spin angular velocity of 
body 1. 

In Appendix [A] we show that, in 5*4, forcing by Y 3 m 
resonates with an inertial mode of the sphere at frequencies 



2 

ui = — 

3 



-m ± 



9 - m' 



1/2 



Q 



(13) 



for —2 < rn < 2. Resonance with the octupolar tide there- 
fore occurs when 



±{0, 0.1700, 0.2981, 0.3922, 0.4444, 0.5099, 
0.8944, 1.1766, 1.3333}. 



(14) 



Similarly, forcing by Y 4 m resonates with an inertial mode at 
frequencies that satisfy the cubic equation 

42 (^) 3+63m (^) 2 - 36 ( 2 - m2 )^ 4m ( n - 2m2 ) = °( 15 ) 

for —3 < m < 3. Resonance with the hexadecapolar tide 
therefore occurs when 

77 

- = ±{0, 0.0970, 0.1770, 0.1920, 0.1940, 0.25, 

0.3273, 0.3540, 0.3840, 0.4550, 0.5, 

0.5580, 0.625, 0.6547, 0.9100, 1.1160, 1.25}. (16) 

It is likely that these resonances persist, albeit with modifi- 
cations, in more realistic models of stars and giant planets. 
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The resonances at smaller values of n/fl may be expected 
to play a role in the tidal synchronization of close binary 
stars and possibly of hot Jupiters, although in the latter case 
tidal potential components beyond the quadrupolar ones are 
intrinsically very weak because of the large ratio of orbital 
semimajor axis to planetary radius. The resonances at larger 
values of n/Q., close to unity, could be important for the more 
rapidly rotating hosts of hot Jupiters, such as F stars, wher e 
spin-orbit misalignments are common l|Albrecht et al.ll2012l) , 
although the convective zones of these stars are of limited 
extent . 



4 TIDAL FORCING OF A ROTATING 
BAROTROPIC FLUID 

4.1 Equilibrium 

A barotropic fluid is one in which the pressure p is uniquely 
related to the density p. This situation arises, in particular, 
when the specific entropy of the fluid does not vary in space 
or time. The differential relations between the pressure, spe- 
cific volume v = 1/p, specific internal energy e and specific 
enthalpy h = e + pv are ther0 de = —p dv and dh = v dp. 
The sound speed is v B = ^Jdp/dp, so that dh — v$ dp/ p. 

In a frame of reference that rotates with uniform angu- 
lar velocity tt, an ideal barotropic fluid satisfies the equation 
of motion, 



~dt 



u ■ Vu + 2SI x u = -V(/i + + $c), 



(17) 



and the equation of mass conservation, 
! + V.(pu) = 0, 

where u is the velocity, <3> g is the gravitational potential, 
which satisfies Poisson's equation, 



(18) 



and 



A-KGp, 



•i|fixr| 2 



(19) 



(20) 



is the centrifugal potential. 

Steady axisymmetric solutions can be found in which 
the fluid is uniformly rotating and therefore has u = in 
the appropriate frame of reference. Such solutions satisfy the 
equilibrium condition 



h + $ g + <I?c = constant. 



(21) 



We assume that the barotropic relation is such that p is 
an increasing function of p for all p > 0, and that p oc p 1+1 / n 
in the limit p — > 0, where n is a (finite) positive real number[f| 
In this case p and p are increasing functions of h for all 
h > 0, and p oc h n , p oc h n+1 and a ft in the limit 
h — >■ 0; the simplest example of such a barotropic relation is 
the polytrope, in which these power laws hold for all h > 0. 



2 If the fluid is barotropic for some other reason, then e and h 
are not the thermodynamic internal energy and enthalpy, but play 
equivalent roles in the fluid dynamics. Some heat exchange with 
the surroundings is then implied. 

3 More precisely, dlnp/dlnp — > 1 + 1/n as p — > 0. 



The body then exists where h > 0, and can have a free 
surface on which h = 0, while the normal enthalpy gradient 
is typically non-zero there. If n < 1, the normal density 
gradient is weakly singular on the surface. 

We will also consider the case of a homogeneous incom- 
pressible fluid, for which p = constant. Although this case 
resembles the limit n — >• of the family of polytropes, it is 
best treated separately because the density is non-zero at 
the free surface. 



4.2 Linearized equations 

We consider a steady axisymmetric body ('body 1') that 
is weakly perturbed by a tidal gravitational potential ^f, 
which need not be steady or axisymmetric. The tidal poten- 
tial is generated by a external body ('body 2') and satisfies 
Laplace's equation V 2 ^ = within body 1. 

The linearized equations governing the response of 
body 1 are 



C + 2fix5 = -VW, 
W = h' + $' + *, 

vV = A-nGp, 



(22) 
(23) 
(24) 
(25) 



where the prime denotes an Eulerian perturbation, the dot 
denotes d/dt, and £ is the displacement, such that u' — £. 
Note that $' represents the internal (self-) gravitational po- 
tential perturbation, while ^ is the external (tidal) poten- 
tial. Outside the fluid, $' satisfies Laplace's equation and 
decays to zero as \r\ — > oo. Equation (|24p implies that 



ti 



(26) 



The desired solution of these equations is such that £, 
V ■ £, W and VW are bounded everywhere within the fluid, 
while $' and V"!?' are bounded and continuous everywhere 
and tend to zero as \r\ — > oo. If n < 1, then p' is weakly 
singular on the surface, but equation (|25[) still has a solution 
such that $' and V"! 1 ' are bounded and continuous. 

The body may contain a solid core. Since a discussion of 
solid mechanics is beyond the scope of this paper, we assume 
that the core is perfectly rigid. The boundary condition £ • 
n — then applies on the surface of the core. However, we 
will also consider a fluid core in Section [49] 



4.3 Love numbers, energy and angular momentum 

From the above linearized equations and boundary condi- 
tions, after some integration by parts, we obtain an energy 
equation of the form 



dE _ 
~dt ~ 
where 



(27) 



is the (canonical) energy of the perturbation, and 



P= / p£- (-V#)dV 



(29) 
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is the power input by the tidal force. Integrals that involve p 
are carried out over the entire fluid volume of body 1, while 
those that do not are carried out over all space. The energy 
E is positive definite for a gravitationally stable body. 
After an integration by parts, we have 



P = 



*V 2 <i>'dV. 



(30) 



The latter integral can be carried out over any region that 
includes the fluid volume of body 1. Let V* be any such 
region that is simply connected and does not include body 2, 
so that V 2 ^/ = in V*. Then, by Green's second identity, 



P = 



1 
1 

4^G 



($'V 2 * - *V 2 $')dV 



($'V* - *V<j>') ■ dS. 



(31) 



av. 



We assume that V* can be chosen such that its bound- 
ary dV* is a sphere of some radius a centred on body 1, i.e. 
that the two bodies are sufficiently distant that they can 
be separated by a spherical boundary. Then, in the vicin- 
ity of dV, , 9 can be represented using an interior multipole 
expansion in solid spherical harmonics, while $' can be rep- 
resented using an exterior multipole expansion. Thus 

* = (£)V(0,0, (32) 



l,m 



Mi) 



-(i+i) 



ir"(M), 



(33) 



where R is the nominal (e.g. equatorial) radius of body 1 
and the sum is over integers I > 2 and —l<m<l. Since *!/ 
and $' are real, we have = (*f )* and = ($")*• 

Having chosen dV* to be the sphere of radius a, we find 



P = 



E 



(2l + l)R 
4ttG 



(34) 



which is independent of the parameter a, as expected. 

An initial-value problem can be solved by Fourier (or 
Laplace) transform methods. If the tidal potential is applied 
for a finite time or decays sufficiently rapidly as t — » oo, 
the total amount of energy transferred to the body is, by 
Parseval's theorem. 



AE 



E 

l,m 

E 

l . m 



where 



Pdt 

(2l + l)R 
4tvG 

(2l + l)R 
AirG 



VT{t) e iut dt, 



*r(*)[*r(*)]*d< 



(-iw)#rH[*r(w)]' ^, 



(35) 



(36) 



etc., denote temporal Fourier transforms. Taking the Fourier 
transform of the linearized equations (|22fl (|25 P leads to a 
purely spatial problem in which the frequency uj (which may 



in general be treated as a complex number) appears as a 
parameter. For a body of arbitrary shape, the solution of 
this problem yields a linear relation of the form 



(37) 



involving an array of potential Love numbers fc™; m (oS). The 
Love numbers are dimensionless complex coefficients that 
quantify the frequency-dependent response of the body to 
tidal forcing; since gravity is the only means of communica- 
tion between the two bodies, all the information relevant to 
spin-orbit coupling and tidal evolution is contained within 
them. If the body is spherically symmetric then fc™; m (cu) — 
unless /' = I and m! = m. Let us assume instead that the 
body is axisymmetric (in the selected spherical polar coor- 
dinate system), so that fc™; m (to) — unless m' = m. Since 
the linearized equations are invariant under complex conju- 
gation, we also have fc ; ™,; m (aj) = [fc ; _ ; 7 l ' _m (— uj*)]* . Then 

A£ = EE (2 ' +1)i? 



4ttG 



j 

(-i^rn^^mrn]*^. 



(38) 



Consider an idealized problem in which only a single 
component ^f 1 of the tidal potential is applied. In fact, 
because * is real, *P;~ m = (^T)* must also be present, 
and their Fourier transforms are related by $^" m (w) = 
[§P(-w*)]*. Then, if m + 0, 



AE = 



+ 



= 2 



(21 + 1)R 
4ttG 

(21 + 1)R 
+ 4ttG 

(21 + 1)R 
4ttG 

(21 + 1)R 
+ 4ttG 

(21 + 1)R 
4ttG 

(21 + 1)R 



(-ioj)krn\*r(co)\ 2 



2tt 



j 

(-MV m HI*r ra HI 2 ^ 



(-iw)fcr(«)i*r(w)r 



271 



c-iw)[*r(-w)]*i*rc-")r ^ 



(-iw)fcT(a;)|*r(w)| : 



do; 
2^T 



4ttG 

(2l + l)R 
AtyG 



(+iw)[*r(w)]i*r(w) 



2tt 



o;Im[fcrH]|*rH| 2 ^, (39) 



where fe™(w) 



! (cj) is the potential Love number that 



quantifies the part of the response that has the same form as 
the applied tidal potential. (If instead m = 0, then the factor 
of 2 is absent from this result.) Similarly, if J represents the 
component of angular momentum parallel to the axis of the 
coordinate system, then the quantity transferred is 



AJ = 2 



(2l + l)R 
4ttG 



mlm[*« m (w)]|*r(w)r^. 



(40) 



If the body is stable, the Love number should be an 
analytic function of u) in the upper complex half-plane. In 
the absence of dissipative processes it may have singularities 
on the real axis, indicating an unbounded resonant response. 
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The Fourier-transformed problem is equivalent to one in 
which the forcing and the response are assumed to depend 
harmonically on time through a common factor e -1 "'. If uj 
is real, and dissipative terms are included in the equations 
to move the singularities below the real axis, then there is 
a steady rate of energy transfer, equal to the rate of energy 
dissipation and proportional to ui\m[k] n (ui)]. Specifically, if 
* = Re[A(r/R) , Y l m (6,</>)e~ iut ], where A is an arbitrary 
constant of appropriate dimensions, then this rate is 

P= {2l ^ R \\A\ 2 ulm[kT^)], (41) 

and the tidal torque on the body differs by a factor of m/uj. 
(In the case m — this expression for P represents a time- 
average.) 

The rate of energy transfer is a frame-dependent quan- 
tity. If the problem is analysed in a frame in which the body 
rotates, then the rate of energy transfer differs from the dis- 
sipation rate. The difference is equal to the rate of change 
of rotational energy of the body under the action of the 
tidal torque. For a similar reason, in a differentially rotating 
body there is no simple global relation between the dissipa- 
tion rate and the rates of energy and angular momentum 
transfer. 

4.4 Low-frequency oscillations in a slowly rotating 
body 

A body of the type considered in Section 14.11 is capable of 
several types of oscillation: acoustic waves, surface gravity 
waves and inertial waves. The frequencies of the first two 
types of oscillation are typically greater than the character- 
istic dynamical frequency (GM/ R 3 ) 1 ^ 2 of the body, while 
the frequencies of the third type of oscillation are typically 
comparable to Q. We define the dimensionless parameter e 
via 




When the body is slowly rotating (e 2 <C 1) and the equi- 
librium structure is close to spherical symmetry, the inertial 
waves are well separated in frequency from the acoustic and 
surface gravity waves. The frequencies of tidal oscillations 
typically lie in the range of inertial waves. 

An asymptotic theory of low-frequency forced os- 
cillations in a slow ly rotating body was developed by 
lOgilvie fc Linl (|2004h . We further develop this theory here 
for the case of barotropic fluids, but without using a formal 
asympt otic expans ion; si milar approximations have been 
used by|Wu| l|2005al lbri and iPapaloizou fc Ivanovl (120051 ). We 
are interested in a low-frequency limit in which the tidal fre- 
quency is comparable to Q, and therefore O(e), with e<l, 
while the dynamical frequenc^Q is O(l). Assuming that the 
tidal potential = O(l) (which is a convenient but arbi- 
trary choice for a linear analysis), we find that the scalings 
£ = O(l), W = 0(e 2 ), h' = O(l), = O(l) are required in 
order to satisfy the linearized equations. To leading order, 
therefore, we wish to solve the system 

4 In effect, we are considering a family of slowly rotating bodies 
labelled by the small parameter e, and adopting a system of units 
such that M and R are the same for each body, while Q oc e. 



i + 2£lx£ = -VW, (43) 

+ + * = 0, (44) 

P = -V • (pO, (45) 

VV = 4ttGp', (46) 



in which the W term has dropped out of equation (|44|) . 

Furthermore, to this level of approximation, the ba- 
sic state can be assumed to be spherically symmetric and 
is conveniently described using spherical polar coordinates 
(r,9,<j>), with the coordinate axis coinciding with the axis 
of rotation. The fluid occupies the region aR < r < R, if 
we assume the solid core to be spherical and of fractional 
radius a with < a < 1, and is in hydrostatic equilibrium 
with dh/Ar = — g, where g(r) is the inward gravitational ac- 
celeration. The boundary conditions are as described above. 
Since the regular solution satisfies b! — g£ r at r — R (see 
equation I26|) . we have 

f r = at r = aR, (47) 
$' + * 

f r = — at r = R. (48) 

g 

Combining equations (|44|l and (146f> , we obtain a type of 
inhomogeneous Helmholtz equation (cf. IZahr]ll966af ). 

VV + ^($' + *) = 0, (49) 

for the potential perturbation within the fluid, which sim- 
plifies to Laplace's equation in the surrounding vacuum and 
in the solid core. For given ^ of an appropriate form, this 
equation has a unique solution such that $' and V<E>' are 
bounded and continuous everywhere and tend to zero as 
\r\ — ► oc. 

In particular, the unforced problem ^ = has the so- 
lution $' = 0, which requires p' = h' = 0, as well as 
£ r = at r = R. Free oscillations in the low-frequency 
limit, which correspond to unforced inertial waves, there- 
fore satisfy the anelastic constraint V ■ (p£) = and the 
rigid boundary conditions £ r = on surfaces with both vac- 
uum and solid. These conditions ensure that the otherwise 
unknown W in equation (|43|l can be determined instanta- 
neously in terms of £, because we obtain a modified Pois- 
son equation V ■ (pVW 7 ) = —V ■ (2pCl x £) with Neumann 
boundary conditions (or regularity at r = R), whose solu- 
tion is unique except for an irrelevant additive constant (or 
function of time). 

In the forced problem $ 7^ 0, we must have $' 7^ 0, 
p' 7^ and h! =fc 0. It is natural to separate the forced so- 
lution into non-wavelike and wavelike parts: £ = £ n w + £w, 
etc. The rationale behind this decomposition is that the non- 
wavelike part is an instantaneous hydrostatic response to the 
tidal potential; this is calculated by neglecting the Coriolis 
force, which is responsible for the existence of inertial waves, 
but it involves the correct Eulerian perturbations p' and 
h! to satisfy the tidal forcing. The wavelike part merely cor- 
rects for the neglect of the Coriolis force. The two parts 
therefore satisfy the following equations, which are a natu- 
ral decomposition of equations (|43[) - (|46fl : 

£nw = -VW nw , (50) 

C + t + * = 0, (51) 
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Pnw = -V • (p£nw), (52) 

V 2 $^ w = 47tG P : iw , (53) 

for the non-wavelike part, together with 

C'w +2flx? w = -VW W + /, (54) 

V ■ (p€w) = 0, (55) 

for the wavelike part with p'„ = /i^ = = 0, where 

/ = -2to x £ nw (56) 

appears as an effective force per unit mass driving the wave- 
like part of the solution. Since f/ v = h' w = $w = 0, we can 
omit the subscript on p' nvi , h' nv , and & nw . 

The boundary conditions are also naturally decomposed 



= V ■ £ c 



£nw,r =0 at r = Q_R, 

Snw.r — 3>t 7* — /t, 

Cw.r = at r = a-R, 

Cw,r = at r = R. 



(57) 

(58) 

(59) 
(60) 



This implies that the wavelike part is driven only by the 
unbalanced Coriolis force acting on the non- wavelike part, 
and not through an inhomogeneous boundary condition. 

Although equations (I50|) (|53[) contain time-derivatives, 
they do in fact imply that £ nw is (or may be assumed to 
be) instantaneously related to the tidal potential ^. This is 
most easily seen by rewriting equation (f50|) as £ nw = -VX, 
where X — W nv ,. The equations and boundary conditions 
defining £ I1W are then instantaneous in time. Note that 



V-(pVX) = p' = -4($' + *). 



(61) 



The right-hand side of this equation is known from the so- 
lution of the Helmholtz-like equation (|49[) . This is, again, 
a well posed elliptic boundary-value problem which deter- 
mines X up to an additive constant. The boundary condi- 
tions (|57[1 and (|58[) correspond to dX/dr = at r = aR and 
regularity of X at r — R. 

This procedure does not provide the most general solu- 
tion of equations (|50|) - ()53[) . because the displacement could 
also contain a vortical part (not satisfying £ = -VX) with 
a linear dependence on time, which would not contribute 
to £. However, any such residual displacement, which would 
depend on the initial conditions, can be considered to be 
part of the wavelike part of the solution. In other words, the 
initial conditions are decomposed such that £ nw is defined 
as above and does not depend on the initial data, while £ w 
contains the two functional degrees of freedom arising from 
the choice of initial values of £ and £. 

The non-wavelike part as defined above is not equi valent 
to the conv entional 'equilibrium tide' as defined by IZahnl 
|l966al ) and lGoldreich fc Nicholson! l|l989l) . The equilibrium 
tide is a spheroidal displacement of the form 



£ c = e r R c — e r x (e r x VS C 



with 



Re = 



(62) 



(63) 



(64) 



It is supposed to represent the low-frequency limit of the 
tidal response, but does not apply in a barotropic body 
where there is no st able stratification, as pointed o ut by 
iTerquem et"afl (|l998l l and lGoodman fc Dickson] (|l998n . The 
equilibrium tide satisfies all equations (|50|) (I53p except the 
first, because it is not irrotational in general. In fact, 



V X £ e = -er X Vh (Re ~ ^) 



and 



rVt R 



dS c 
dr 



= V 2 (rR e ) = -V 2 



-($' + *) 



(65) 



(66) 



which does vanish in a homogeneous body, in which jar 
and V 5>' = 0, but not generally. In addition, if the body 
contains a (perfectly rigid) solid core, the equilibrium tide 
does not satisfy the appropriate rigid boundary condition 
there. 

The equations governing the wavelike part of the so- 
lution are well posed as an initial-value problem, even in 
the absence of dissipative processes. In contrast, when the 
forcing is harmonic in time and the solution is assumed to 
have the same periodicity, the equations governing its spatial 
structure are hyperbolic and not well posed in the absence 
of dissipation for frequencies in the range — 2Q, < oj < 2f2 of 
inertial waves, and singularities may occur. 

The energy equation for the wavelike part is 



(67) 



where the integrals are over the fluid volume. The right- 
hand side represents the power input to the wavelike part 
from the effective force; the energy comes ultimately from 
the source of the external potential. Only kinetic energy is 
relevant here because of the absence of stable stratification 
and the inability of inertial waves to elevate the free surface 
in the low- frequency limit. If we introduce viscous or fric- 
tional forces to provide a mechanism of dissipation and to 
resolve the singularities in the response, then a dissipative 
term should be included in this energy equation. 



4.5 Fourier transform and harmonic forcing 

The initial-value problem for the wavelike part can be anal- 
ysed by taking a Fourier (or Laplace) transform in time, 
leading us to consider the problem 



(68) 
(69) 

(70) 



etc. The rigid boundary conditions £ w , r = apply and we 
may wish to consider complex frequencies uj in order to carry 
out the inverse transform and to avoid singularities. The 
Fourier-transformed variables u„, etc., are also complex in 
general. 



— ium w + 2f2 x u w = — VW W + /, 

V • (pu w ) = 0, 

where u w = — kj£ w and 

£ w (r,w)= / ^ w (T',t)e 1 "*dt, 
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If the force is applied for a finite time or decays suffi- 
ciently rapidly as t — > oo, the total amount of energy trans- 
ferred to the wavelike part is, by Parseval's theorem, 



AE = 



pu„ ■ f dV dt 



pu w ■ / dV 



2tt' 



(71) 



The Fourier-transformed problem is formally equivalent 
to one in which the forcing and the response are assumed 
to depend harmonically on time through a common factor 
e _ , with lu being complex in general. In this case the 
interpretation of the equations is different, however. The 
force and the wavelike displacement are Re[/(r) e -1 "'] and 
Re[£ w (r) e _iwt ], respectively, and (when lv is real) the aver- 
age rate of energy transfer is 



P = - Re 

2 



pu„ ■ / dV". 



(72) 



As noted above, this problem is generally ill posed for real 
frequencies in the range of inertial waves, so we may choose 
to add a viscous or frictional force to ensure a regular solu- 
tion. In that case the dissipation rate equals the power input 
for real frequencies. For example, if the frictional term — 7ti w 
is added to the right-hand side of equation (|68[) , where 7 > 
is a frictional damping rate, then we obtain 



[Im(u) + 7] / p|u w | 2 dV = i Re / pu w ■ /* dV. 



(73) 



We quantify the solution of the harmonic forcing prob- 
lem as follows. If the force derives from a tidal potential 
* = Re[A{r/R) l Y l m {9,<j})e- iujt ], where A is an arbitrary 
constant of appropriate dimensions, we let 



p£w/ dV = 



(2l + l)R lAl2 
4ttG lA 



(74) 



thereby defining the complex dimensionless response func- 
tion if™ (a;), which is linearly related to the response 
of the fluid and has the Hermitian symmetry property 
KT(u) = [K- m {-u)*)]* . This is defined by analogy with 
equation (|4ip , which relates the energy transfer to the 
potential Love number fc ; m (t<j). Comparing equations (|72[) 
and (I74|l . we find that, for real frequencies, 



= (2l + l)R l 2 
4nG 2 1 1 



(75) 



which agrees with equation (|4Tj) if Im[fc" l (^)] = {to)]. 
Indeed, in the low-frequency limit considered here, we ex- 
pect that fc ; m (cj) is real to a first approximation, and that 
the dissipative part of the response is given by Im[fcf (w)] = 
Im[A"; m (cj)] = 0(e 2 ). Since the effective force driving the 
wavelike tide is a Coriolis force, / oc Q, the response £ w is 
also proportional to fl, and therefore K™ oc Q, 2 . The reason 
for quantifying the response using equation Q74p rather than 
equation (|37|) is that, in the low-frequency limit, the gravi- 
tational potential perturbation $' is a quantity that is not 
readily available, hidden as it is within W . 

When a single potential component \I f ] 71 (t) (having a 
general time-dependence that decays sufficiently at large t, 
and being accompanied by ^ ; -m when m ^ 0) is applied, the 
total energy transferred to the wavelike tide is, by analogy 
with equation (|39[) . 



A£ = 2 ^4^ r-Im[A7H]|*rH| 2 ^, (76) 
J — 00 

or half this quantity in the case m = 0. 



4.6 Impulsive forcing 

We now consider a special tidal forcing problem in which 
the wavelike part experiences an impulsive effective force of 
the form 



f = f(r)S(t), 



(77) 



where S(t) is the Dirac delta function. Since / is related 
to £nw through equation (|56[l and £ n w is instantaneously 
related to the tidal potential, this means that the potential 
should be of the form 



(78) 



where H(t) is the Heaviside step function. This, rather than 
the more obvious $ oc 5(t), is the appropriate type of im- 
pulse problem to consider within the low-frequency approx- 
imation. The impulse is supposed to occur slowly enough 
that the anelastic constraint holds, i.e. slowly compared 
to the sound crossing time, but fast compared to the ro- 
tation period. The idea is to excite a broad spectrum of 
inertial waves but no surf ace gravity or acoustic waves. 
IPapaloizou fc Ivanovl (|2010h in their Appendix BI discuss a 
problem of this type, but without explicit calculations; they 
do, however, solve initial-value problems in which inertial 
waves are excited by a parabolic tidal encounter. 

Assuming that the fluid is at rest before the impulse, 
immediately afterwards it will have a wavelike velocity given 
by 



Uy, = f - VWw, 



(79) 



where W v is arranged to satisfy the anelastic constraint 



V ■ (pu w ) = 



(80) 



and the boundary conditions u„, r = 0. Note that iiw is 
just the initial velocity immediately after the impulse; it will 
subsequently oscillate as a collection of inertial waves and be 
damped, if dissipative terms are included in the equations. 
Equation (|79[) is derived by integrating equation (|54[) over 
an arbitrarily small time-interval that includes the instant 
t — 0. Note that the Coriolis force acting on the wavelike 
velocity has no effect during the impulse process; the only 
role of the Coriolis force is to provide the effective force / 
from its action on the non-wavelike velocity. 

We now require a procedure to calculate the impulsive 
response and the associated energy transfer 

PKfdK (81) 

We consider a tidal potential with the spatial structure 



* = Re *((r)Y," 



Re 



A 



(£)V<*. 



(82) 



The associated internal gravitational potential perturbation 
[also proportional to H(t)], 



$' = Re [M(r)Yr(e, 



(83) 
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satisfies the Helmholtz-like equation (|49|) . which reduces to 
the ordinary differential equation (ODE) 



J__d 

r 2 dr 



$; + l^ ( $; + * i) = (84) 



in aR < r < R. Matching to decaying solutions of Laplace's 
equation in the solid core and the exterior vacuum provides 
the boundary conditions d<J?(/dr = l&i/r at r — aR and 
d$J/dr = —(I + l)$j/r at r = R. The associated non- 
wavelike tide [also proportional to H(t)] is 



•VI, 



(85) 



where X = Ke[Xi(r)Y l m (9, (j>)\ satisfies equation ()6ip . which 
reduces to the ODE 



_l__d / 2 dX, 

r 2 dr \ ^ dr 



1(1 + 1) 



pX t 



(86) 



in aR < r < R, together with dX;/dr = at r — aR and 
regularity of Xi at r — R. From £ nw we deduce the impulsive 
effective force / = — 2tt x £ nw , which has both spheroidal 
and toroidal parts. Using standard methods of pro jection on 
to vector spherical harmonics (e.g. lOgilvie fc Linll2004h we 
find from equation (|79[) that the impulsive wavelike velocity 
is of the form 

u w = Re [e r a ; (r)y ; m + r 2 b,(r)VYi m 

-r 2 Q_i(r)e r x V*Ti - r 2 c l+1 (r)e r x VY^] , (87) 

where &i and b\ represent the spheroidal part and &i±\ the 
toroidal part. These components satisfy 



o; = 



h = 



CJ-1 



Ci + l 



2imQ.Xi 


_ dWi 


r 


dr ' 


2imSl 


- (r dJCl 


l(l + l)r 


2 \ dr 


2Qqi 


r— — + ( 




dr 




( dX t 







+ x l] -^ 



(I + l)Xi 



IX, 



\^-(r 2 pai)-l(l + l)pbi = 0, 
r A dr 

where W v = Re[W' ; (r)F ! m (6», <t>)} and 

1/2 



<ln 



1 / n z -m z 
n \ 4n 2 - 1 



(89) 

(90) 

(91) 
(92) 

(93) 



is a coefficient arising from the coupling of spheroidal and 
toroidal velocity components by the Coriolis force. Therefore 
Wi satisfies 



l__d / 2 dWi 

r 2 dr \ ^ dr 



1(1 + 1) 



pWi = 



2imf2 dp ~ 

1~ Xl > 

r dr 



(94) 



with boundary conditions dWi/dr = — 2\m£lXijr (corre- 
sponding to a, = 0) at r = aR and r — R. 

The impulsive energy transfer is then found to be 



E = Ei + Ei-i + Ei+i, 
with 



(95) 



Ei = \ / pr 2 [\ai\ 2 + l(l + l)r%\ 2 ] dr, 

J aR 

i r R 

Ei-i = - / pr 2 [i(I-l)r 2 |ft_i| a ] dr, 

J aR 

i r R 

El+l = j / pr 2 [{l + l)(l + 2)r 2 \ci +1 \ 2 ] dr, 

J aR 



(96) 
(97) 
(98) 



although these expressions should be doubled in the case 
m — 0, assuming that A is real. Let us now compare this 
result with equation ([76[1 , The impulsive (switched-on) po- 
tential corresponds tq_| 



1 y ' 2uj' 
and so 

&=2 (2i+i)r\ai r Im[i , rM] ^_ 

AttG Sty 1 K n u 



(99) 



(100) 



4.7 Interpretation 

The previous subsection shows that the energy transferred 
to the wavelike tide in the impulse problem considered above 
corresponds to a certain frequency- average of the imaginary 
part of the Love number of the body. Although the integral 
as written is over all values of u, in fact we need integrate 
only over the range —2Q. < lo < 2Q where inertial waves oc- 
cur, and where the low-frequency approximation holds0 If 
no dissipation is included in the equations, if; m (w) has singu- 
larities on the real uj axis and the integral can be computed 
only by deforming the integration contour around them in 
an appropriate way. In that case the integral describes the 
energy given to waves that are never damped. If, however, 
some dissipation is included, then the waves are ultimately 
damped and the integral describes the energy dissipated as 
well as that transferred. 

We see, therefore, that the dimensionless quantity 
Im[_K) m (aj)] duj/co, which is a frequency-averaged mea- 
sure of the dissipative properties of the body in the low- 
frequency range corresponding to inertial waves, can be di- 
rectly related to the energy transferred in an impulsive in- 
teraction. Furthermore, this quantity can be computed di- 
rectly by solving a simple set of ordinary differential equa- 
tions. This procedure is much easier than computing an iner- 
tial wave with a harmonic time-dependence, because during 
the impulsive interaction the Coriolis force has no time to 
act on the wavelike velocity; the only role of the Coriolis 
force is to provide the effective force from its action on the 
non-wavelike velocity. The solution depends on the inter- 
nal structure of the body but not on its dissipative prop- 
erties; this solution is smooth and free from boundary lay- 
ers. Following the impulse, this velocity field will resolve 



5 It may be better to think of the impulse as H (t) e et with 
e \ 0. Then *|"(w) = L4/[2(w + ie)], again with e \ 0. 

6 Note that, if a low-frequency approximation is not made, and an 
impulsive gravitational force arising from a tidal potential oc 8(t) 
is applied, then the energy transfer is instead proportional to 
J Im[fc ; m (oj)] udu. This integral is expected to be dominated by 
surface-gravity modes of higher frequency. 
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into a collection of inertial waves and be damped, if dis- 
sipative terms are included in the equations, but for our 
present purposes we need not consider this subsequent evo- 
lution. In the following subsections we compute the quan- 
tity /_ lm[K{ n (uj)] (Iuj/lj for some simple interior mod- 
els and compare with numerical solutions of the frequency- 
dependent response functions of such bodies. 

4.8 Application to a homogeneous fluid 

We first consider, as in lOgilvid l|2009h . a body consisting of 
a homogeneous incompressible fluid of density p, possibly 
containing a perfectly rigid solid core of fractional radius a. 
The mean density of the body is p. 

In this case the Helmholtz-like equation (|84l) becomes 



1 d 



,d$; 



r 2 dr \ dr 
4nGp 



1(1 - 



(101) 



where g = GM/R 2 = 4TvGpR/3 is the surface gravity. Note 
that, in this incompressible model, the sound speed is infi- 
nite, and the Eulerian density perturbation is given by 

p =-(i-Vp = irp8{r-R). (102) 

Given = A(r/R) 1 , the solution is ${ = B(r/R) 1 for < 
r < R and $J = B(R/r) l+1 for r > R, where (21 + l)B = 
(3p/p)(A + B). 

Since V ■ £ = for an incompressible fluid, X satisfies 
Laplace's equation and the solution satisfying the boundary 
conditions ((57} and (1581) is 

(i)' + (nhK +, (ffT « 



X l =C 

with 

C 



(A + B)R 



1(1 - aP'+^g' 



(104) 



Ww also satisfies Laplace's equation and the relevant 
solution is 



Wi = 



2imfiC 
I 



© -Gil) 



Thus 
ai — bi — 0, 
2Qq t 



r 

2Q.q. 



eta + D(;)', 
t^ + D^) 



21 + 1 



Cl-1 = 

h+i = - 

and so 
E t = 0, 

Ei^ = ptt 2 q 2 \C\ 2 l(l - l)(2l + 1)(1 - a 2l+1 )R, 
-.2 ~2 i/-,i2ia ( <• + 2 



r . 



(105) 

(106) 
(107) 

(108) 

(109) 
(110) 



Ei+i = pQ. q~i+i\C\ I 



l + l 



(2l+l)a 2l+1 (l-a 2l+1 )R.(Ul) 



The total impulse energy E is the sum of these three ex- 
pressions. Now qf oc (I 2 — m 2 ) vanishes for sectoral harmon- 
ics (|m| = I), while qf, 1 is non-zero. For sectoral harmon- 
ics, therefore, the impulse energy has a strong dependence 



on core size, being proportional to a 2l+1 for small a. For 
tesseral harmonics (\m\ < I), however, E remains significant 
in the limit a — ► because of the contribution of Ei-\. We 
also see that the impulse energy is proportional to SI 2 , as 
expected. 

For I — m — 2 we have 



E 



^(l-a 5 )\C\ 2 P Rn 2 



20 5 
= 189° (1 



^\A + B\ 



, P R 3 n 2 



(112) 



which is proportional to a° for small a. This 
strong dependence on core size agrees with the 

behaviour of t he (frequency-avera ged) dis sipation 

rate found by iGoodman fc Lacknerl (120091 ') . lOgilvid 
J2009I) and iRieutord fc Valdettard l|2010h . although 
I Goodman fc Lacknerl (|2009l ) considered a nonlinear mecha- 
nism involving wave breaking. In the simplest case p = p, 
we have (A + B)/A= (21 + 1)/(2Z - 2) = 5/2 and so 



t r&-2/ m dw 1007T 2 

Im \K 9 (uj)\ — = e 

1 2V uj 63 



(113) 



We have verified this result numerically as follows. In Fig- 
ure[T]we show the frequency-dependent response to I — m = 
2 tidal forcing for a homogeneous body with an incompress- 
ible fluid envelope and a perfectly rigid solid core of v ariou s 
sizes. The tidal response is computed as in lOgilvid (2009) 
and this figure is equivalent to (parts of) Figures 1 and 2 
of the earlier paper except that the energy dissipation rate 
is converted into the imaginary part of the Love number. 
This conversion brings in the parameter e, defined in equa- 
tion (|42|) . which is unspecified in these calculations but is 
assumed to be small; as noted above, we expect Im fc oc e 2 
for inertial waves. Then the integral in equation (|113|) is 
computed numerically after subtracting the baseline dissipa- 
tion rate for each curve (which corresponds to the frictional 
damping of the non- wavelike part of the tide). The results 
are plotted in Figure [2] for nine values of a and three values 
of the frictional damping coefficient. Very good agreement 
is found with equation (|113|) , 

The divergence that occurs as a — > 1 deserves some 
comment. The coefficient C diverges in this limit, as do 
the non-radial components of the non-wavelike displace- 
ment. The thin, incompressible fluid shell is squeezed radi- 
ally by the tidal force and produces rapid horizontal motion. 
Through the Coriolis force, a large wavelike velocity is also 
generated. This behaviour would be weakened if the core 
were not perfectly rigid but instead underwent some tidal 
deformation. 

Generalizations of equation (|113p can be found in Ap- 
pendix [B] In particular, 



Im[K 2 (u)] — = — e 
cj 504 



and 



Im[xS(w)] 



did 



5n 2 
— e 
14 



189 + 256q j 
1 - a 5 



7 + 8a 5 



(114) 



(115) 



Note that, although these expressions are increasing func- 
tions of a, they do not vanish at a = 0. We have also 
found excellent agreement between these expressions and 
the numerically integrated Love numbers using the method 
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Figure 1. Frequency-dependent tidal response to Y 2 2 of a homogeneous body with an incomp ressible fluid e nvelope and a perfectly 
rigid solid core of fractional radius a = 0.1 (top panel), 0.3, 0.5, 0.7 or 0.9 (bottom panel). As in lQg ilvie (200^) the fluid is inviscid but 
experiences a frictional force with a damping coefficient 7 given by 7/n = 10 — 2 (red line), 10 — 2 6 (green line) or 10~ 3 (blue line); also, 
the surface is not free but has a prescribed radial motion. The imaginary part of the Love number is deduced from the dissipation rate 
and is plotted on a logarithmic scale after multiplication by e -2 (fi / 'a) , where f! = t(GM / ' R?) 1 / 2 . A numerical resolution of L = N = 400 
is used. 
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Figure 2. Frequency-integral of the imaginary part of the Love 
number, weighted by 1/u and divided by e 2 , for a homogeneous 
body, versus fractional core size. The solid line shows the ana- 
lytical result II113H of the impulse calculation, while the symbols 
show the numerical integrals based on Figure [l] and similar cal- 
culations. 



described above for the case I = m = 2. The frequency- 
dependent responses are illustrated in Figures [3] and [4] Note 
that the enhancement of the tidal response in these cases 
results from resonances with the spin-over mode (m = 1, 
uj/Q. = —1) and the spin-up mode (m = 0, ui/fl = 0). These 
are trivial modes and there should be no associated dissi- 
pation. In these calculations the dissipation is artificial and 
results from the use of a frictional, rather than viscous, force. 
In either case, however, there is a formal singularity in Imfe 
that gives a non- vanishing contribution to the integral even 
when a = 0. In a viscous fluid the resonance should be of 
zero width. 

To illustrate the richer possibilities that exist for tidal 
potential components beyond the quadrupolar ones, we show 
similar results for I = 4 and m = 1 in Figure [S] In this case 
the response for small values of a is dominated by reso- 
nances with non-trivial inertial modes of the full sphere (at 
w/n = -1.7080, -0.6120 and +0.8200, as expected from 
equation ll5[) . For larger core sizes the frequency-dependence 
is qualitatively similar to that seen for quadrupolar tidal po- 
tentials. The analytical result in this case is 



Tm[Kl(u)] — 



27n 2 f 6875 + 7168a s 



616000 V 



1 



(116) 



The integrated responses for all tidal potentials up to 
I — 4 are shown in Figure [6] Although there is a general 
weakening of the response with increasing I, the integrated 
response to a tesseral harmonic of higher I is usually much 
stronger than that to a sectoral harmonic of lower I. 



4.9 Application to a piecewise-homogeneous fluid 
model 

We now examine the effect of replacing the rigid solid core 
with a fluid core. We therefore consider a homogeneous fluid 



of density p\ in < r < aR and another of density p 2 in 
aR < r < R, with pi > p2 for gravitational stability. 

At an interface, where the density changes abruptly, 
£nw,r = — ($' + 9)/g. This is because the Lagrangian pres- 
sure perturbation p' — pg£ r is continuous, but p' = p(W — 
$' — \E') and W is supposed to be negligible [0(e 2 )] in the 
low-frequency approximation. If p is discontinuous then we 
find the stated condition, which means that the interface 
moves equipotentially. (This should be valid provided that 
the tidal frequency is small compared to the relevant inter- 
facial gravity wave frequency.) 

The Helmholtz-like equation is then 

V 2 $'+^($'+*)[(p 1 -p 2 )5(r-a7?)+p 2( S(r-J?)] = 0.(117) 
Given * = AirjKfY^, the solution is $' = $j(r)Y; m , with 



(B 1+ B 2 )(r/R) L , 

B^r/R) 1 + B 2 a 2l+1 (R/r) l +\ 

(B 1 +B 2 a 2l + 1 )(R/r) l+1 , 



< r < aR, 
aR < r < R, 
r>R, 



(118) 



where the coefficients satisfy the matching conditions 

(2l + l)B 2 =3(1- f)(A + B 1 +B 2 ), (119) 



(21 + l)Bi 



3/ 



T (A + B 1 +B 2 a 



2i + l\ 



(120) 



f + (i-fW 

with / = p 2 /pi- We shall not write explicitly the algebraic 
solution for Bi and B 2 . 

X satisfies Laplace's equation in each fluid shell and the 
solution (regular at r = 0) is of the form 



X, = 



Cx(r/R)\ 0<r<aR, 
C 2 (r/R) 1 + C 3 a a+1 (R/r) l+1 , aR < r < R. 



(121) 



The condition on the radial displacement is dXi/dr — ($; + 
^i)/g at each interface. Thus 



Id = IC 2 - (I + 1)C 3 
IC 2 -(1 + l)C 3 a 2l+1 = 



_ 3(A + 5i + B 2 ) 
AnGp! 

3(A + gi + B 2 a 2l+1 ) 
4:nG(p 2 + (pi - p 2 )a 3 ) ' 



(122) 



(123) 



The fact that C3 7^ in general means that Xi is discontinu- 
ous in general at r — aR. This means that the non-wavelike 
displacement has a tangential discontinuity. The singularity 
of vorticity results from the non-barotropic nature of the 
density jump. 

W also satisfies Laplace's equation and the solution is 



Wi -- 

such that 
&i = bi = 0. 



-2imQ(Ci/l)(r/R) 1 , 

-2imfl[(C 2 /l)(r/R) 1 - (C 3 /(l + l))a 2l+1 (R/r) l+1 ] 



(125) 



(Here and below, the first and second cases of the expres- 
sions in braces apply for < r < aR and aR < r < R, 
respectively.) The toroidal part is given by 



C;-i = 



C1+1 — 



-2nq l r- 2 (2l + l)C 1 (r/R) 1 , 
-2nq l r- 2 (2l + l)C 2 (r/R)\ 



0, 



(126) 
(127) 



-2flqi +1 r- 2 (2l + l)C 3 a 2l+1 (R/r) 
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(124) 




The contributions to the impulse energy are then 

E x = 0, (128) 

Ei-i = l(l~l)(2l + l)Q 2 qf 

x[ Pl a 2i+1 |d| 2 + P2 (1 - a 2l+1 )\C 2 \ 2 ]R, (129) 



E l+1 = (l + l)(l + 2)(2l + l)Q 2 q 2 +1 

xp 2 a 2l+1 (l-a 2l+1 )\Cs\ 2 R. (130) 

These results are converted into integrals of Im k using 
equation (|76[) and plotted in Figure [7] for various values of / 
in the case I — m — 2; the analytical result is given in 
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Figure 4. As for Figure [T] but for the tidal potential ■ There is a formal resonance with the spin- up mode at ui/Q, = 0. 



equation (|B3|) . For small values of / (large density contrasts) 
the impulsive response is similar to, but weaker than, the 
case of a homogeneous fluid with a solid core; the differences 
are greater when the core is larger. These results show that a 
fluid core can play a similar (if weaker) role in the excitation 
of inertial waves to that of a rigid solid core. As / tends to 1 



the response disappears because the system becomes a full 
homogeneous sphere. 
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4.10 Application to polytropes 

We have also computed the impulsive energy transfer E for 
polytropes of various indices n. This is done by first solving 
the Lane-Emden equation for the structure of the polytrope, 
then computing the impulsive response by solving the ODEs 



of Section 14.61 numerically by standard methods. We allow 
for a perfectly rigid solid core to be present, and assume 
that the mass of the core is equal to the mass of fluid that 
would occupy the equivalent region in a full polytrope. 

The frequency-dependent response of an n = 1 poly- 
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Figure 8. Frequency-dependent tidal response to of an n = 1 polytrope with a fluid envelope and a perfectly rigid solid core of 
fractional radius a = 0.1 (top panel), 0.3, 0.5, 0.7 or 0.9 (bottom panel). The fluid is viscous and the Ekman number is Ek = 10 -3 (red 
line), 10 — 4 (green line), 10 -5 (blue line) or 10 — 6 (black line). The imaginary part of the Love number is computed directly and is plotted 
on a logarithmic scale after multiplication by e~ 2 (Q/uj), where ft = e(GM/R s ) 1 ' 2 . In these calculations e = 0.01 and the fluid envelope 
is truncated with a free surface at a fractional radius of 0.99. A numerical resolution of L = N = 200 is used. The results for Ek = 10 — 6 
are not fully converged at this resolution and occasionally return negative values of e _2 (f2/oj)Im k, resulting in gaps in the black line. 
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Figure 6. As for Figure [2] but showing only the results of the 
impulse calculation for a homogeneous body and for various tidal 
potentials with I = 2 (red lines), 1 = 3 (green lines) and I = 4 
(blue lines). On the left of the figure the ordering of curves is, 
from top to bottom: Y° Y} , Y° Y,\ Y 2 Y° Y.\ Y? , Y? , 
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Figure 7. As for Figure [2] but showing only the results of the 
impulse calculation for a homogeneous body (black line) and for 
the piecewise-homogeneous fluid model with density ratio / = 
0.25 (red line), / = 0.5 (green line) or / = 0.75 (blue line). 



trope to I = m = 2 tidal forcing is illustrate d in Figure El 
These results are similar to those presented bv lOgilvie Sz Lin! 
l|2004l ) but are computed by a different numerical method 
which will be described in forthcoming work (Ogilvie, in 
preparation). Rather than making a low- frequency approx- 
imation as in lOgilvie fc Linl (|2004h . we solve the full set of 
linearized equations for a compressible, self-gravitating and 
viscous fluid while neglecting centrifugal effects (the parame- 
ter e is here set to 0.01). This method has the advantage that 
the Love number can be measured directly, because the grav- 
itational potential perturbation is part of the numerical so- 
lution. The fluid is subject to viscous, rather than frictional, 




Figure 9. As for Figure [2] but for an n = 1 polytrope. The 
numerical integrals are based on Figurc|8]and similar calculations. 




Figure 10. As for Figure [2] but showing only the results of the 
impulse calculation for a homogeneous body (black line) and for 
polytropes of indices n = 1 (red line), n = 3/2 (green line) and 
n = 3 polytrope (blue line) . 



forces, as quantified by the Ekman number Ek = u/(2QR 2 ), 
where v is the kinematic viscosity. In all cases except those 
at the lowest Ekman number Ek = 10 -6 , where the numeri- 
cal resolution may be insufficient, the imaginary part of the 
Love number is consistent with the viscous dissipation rate 
to a high degree of accuracy. 

While there are many qualitative features in common 
between Figures [T] and [H] there are also important differ- 
ences. For very small core sizes such as a — 0.1, the ho- 
mogeneous body exhibits almost no response to Y 2 , while 
the n = 1 polytrope shows a few resonant peaks. These 
correspond to the resonant excitation of inertial normal 
modes of the full polytrope, which have been computed by 
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lLockitch fc Friedmanl fl999)Q Similar modes also exist in 
the homogeneous sphere, but they have no overlap with 
the y 2 2 tidal potential and are not excited. For larger core 
sizes, normal modes cease to be important, until the limit 
of a thin shell is reached, in which the dominant feature is 
the excitation of a planeta ry wave/Rossby wave/r mode at 
w/fi = -1/3 (|Ogilvidl2009l l. 

In Figure [9] we show the integrated quantity 
J Im k dui/uj, divided by e 2 , as determined both from numer- 
ical integration of the curves in Figure [8] (and similar calcu- 
lations) and from the impulse calculation described above. 
The two are in very good agreement, except in a few cases 
where the numerical integration probably has insufficient 
frequency-resolution to capture the narrow peaks that oc- 
cur in Imfe at low Ek, or where the numerical resolution is 
insufficient to determine the dissipation rate accurately. 

Figure [10] compares similar results (but without the nu- 
merically integrated values) for a homogeneous body and 
for polytropes of indices n = 1, 3/2 and 3. This figure shows 
that, as the polytropic index increases and the body becomes 
more centrally condensed, the integrated tidal response to 
y 2 2 becomes less sensitive to the size of the core. The main 
reason for this result is that the overlap between the Y% po- 
tential and the inertial normal modes of the full polytrope 
increases as n increases, so that a significant tidal response 
(albeit concentrated into narrow peaks) is possible even in 
the absence of a core. 

These results are in broad agreement with the findings 
of IPapaloizou fc Ivanovl (|201[f ). who calculated the energy 
transferred (predominantly to inertial waves) to a slowly ro- 
tating n — 1 polytrope during a distant parabolic encounter. 
They found that the energy transferred via Y% was increased 
by nearly an order of magnitude when a (perfectly rigid) 
solid core with a = 0.5 was introduced, but hardly at all 
when a = 0.25. This is broadly in line with the red curve in 
Figure QJj] although a parabolic encounter provides a tidal 
force that is localized in time, their problem is not identical 
to the impulse problem considered in this paper. For exam- 
ple, the Fourier transform of the Heaviside function has a 
singularity at zero frequency that is not present in the case 
of a parabolic encounter. 



5 CONCLUSIONS 

In this paper we have examined the linear response of slowly 
rotating, neutrally stratified astrophysical bodies to low- 
frequency tidal forcing. This problem may be regarded as 
a simplified model of tides in convective regions of stars 
and giant planets. The response can be separated into non- 
wavelike and wavelike parts, where the former is related in- 
stantaneously to the tidal potential and the latter may in- 
volve resonances or other singularities. The imaginary part 
of the potential Love number of the body, which is directly 
related to the rates of energy and angular momentum ex- 
change in the tidal interaction and to the rate of dissipa- 
tion of energy, may have a complicated dependence on the 
tidal frequency. However, a certain frequency-average of this 

7 Visible in the top panel of Figure [8] for example, are the fre- 
quencies u/Q = 0.5566 and —1.1000 listed in the fifth column of 
their Table 6, albeit with a different sign convention. 



quantity is independent of the dissipative properties of the 
fluid and can be determined by means of an impulse calcula- 
tion. The result is a strongly increasing function of the size 
of the core when the tidal potential is a sectoral harmonic 
(|m| = I), such as the I — m = 2 potential usually con- 
sidered, especially when the body is not strongly centrally 
condensed. However, the same is not true for tesseral har- 
monics (|mj < I), which receive a richer response and may 
therefore be important in determining tidal evolution even 
though they are usually subdominant in the expansion of the 
tidal potential. We have also discussed analytically the low- 
frequency response of a slowly rotating homogeneous fluid 
body to tidal potentials proportional to spherical harmon- 
ics of degrees less than five. Tesseral harmonics of degrees 
greater than two, such as are present in the case of a spin- 
orbit misalignment, can resonate with inertial modes of the 
full sphere, leading to an enhanced tidal interaction. Similar 
behaviour can be expected in more realistic models. 

The calculations carried out in this paper are based on 
linear theory and employ highly simplified interior structure 
models. There are various reasons why the complicated de- 
tails of the frequency-dependent response functions might 
not be applicable to astrophysical bodies. The propagation 
of inertial waves may be impeded by nonlinearity (leading to 
instability, wave breaking, etc.), or interrupted by interac- 
tions with convection, magnetic fields, differential rotation, 
buoyancy forces, etc. The waves may not reflect perfectly 
from the boundaries of the convective regions that support 
them. These complications are likely to tend to wash out 
some of the structure in the frequency-dependent response 
curves. Nevertheless, the integrated response functions cal- 
culated in this paper are likely to be much more robust; 
they can be related, as we have seen, to the energy trans- 
ferred in an impulsive interaction, which is independent of 
the detailed way in which the waves subsequently propagate, 
reflect and dissipate. The normal-mode resonances that we 
have identified with higher-order tidal potentials are also 
likely to be important even in the presence of effects such as 
convection. 

Tidal interactions at low orbital eccentricity can be 
studied by analysing the tidal potential into a small number 
of components with a harmonic dependence on time. In this 
case the integrated responses discussed in this paper can be 
taken as indicative of the typical level of tidal dissipation, 
neglecting the complicated frequency-dependence associated 
with the propagation of inertial waves in a spherical shell. Of 
particular importance are the scaling of the imaginary part 
of the Love number (or the reciprocal of the tidal quality 
factor) with the square of the dimensionless rotation rate e, 
and its dependence on the size of a solid or fluid core that 
excludes the inertial waves by total or partial reflection. At 
high orbital eccentricity, tidal interactions are more impul- 
sive in character because the tidal force is strongly peaked 
near pericentre, and the integrated responses are more di- 
rectly applicable to such situations. 

Previous works investigating tidally forced inertial 
waves in astrophysical bodies have tended to emphasize ei- 
ther global normal modes or the singular phenomena as- 
sociated with wave attractors and critical latitudes. This 
paper attempts to bridge the gap between these descrip- 
tions. We consider that inviscid inertial waves in spherical 
(or spheroidal) geometry form true normal modes only un- 
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der special conditions, i.e. when a core is absent and when 
the density profile is sufficiently smooth. We find classical 
resonances with normal modes in coreless homogeneous bod- 
ies and polytropes. As a core is introduced and increased in 
size, some 'memory' of these modes is retained initially, but 
later the singular wave phenomena dominate. 

The application of these findings to extrasolar planets 
and solar-system bodies requires further work. The interior 
structure of giant planets is still quite uncertain; there may 
be non-smooth features such as discontinuities in the density 
or its radial derivative. Progress in the modelling of plane- 
tary interiors will certainly assist in the determination of 
their tidal responses. 
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APPENDIX A: RESPONSES OF A HOMOGENEOUS SPHERE 



In this section we give some explicit solutions for the linear response o f a slow l y rota ting homogeneous fluid body without 
a solid core to low-frequency tidal forcing. The notation follows that of Ogilvie (2009), with a time-dependence e _mjt in the 
frame rotating with the body, but viscous and frictional damping forces are omitted, so uj n = u + 2mQ/[n(n + 1)]. 

Given * = A(r/R) l Yi n , the solution of the Helmholtz-like equation is, as in SectiongTHl "I"' = B^/RfY™ for r < R and 
$' = B(R/r) l+1 Y{ n for r > R, with 2(1 - l)B = 3A. Then X = C(r/R) l Y, m with C = (A + B)R/(lg) and g = GM/R 2 . This 
gives £nw,r = aR(r / R) l ~ 1 Y{ n , where a = —IC/R 2 is the dimensionless non-wavelike radial tidal amplitude at the surface. The 
non-wavelike velocity is given by 



a, 



-iwaR l 



\r) • 6!= - lw ibJ ■ 

The wavelike velocity solutions in various cases are as follows. 



(Al) 



Al The case I = 2 

5aq 2 tvfl 



Cl 



W 



5aqiq 2 ojQ. 2 
r , 



au>fi(mcji — 5q 2 fl) 2 



(A2) 



2oji wi 2oji 

Note that a constant ci, as occurs here, implies a special type of solution, not really an inertial wave. If m = 2 this term 
vanishes because then (72 = 0. If m = 1 the term ci is a spin-over mode (tilting the rotation axis). If m = it is a spin- up 
mode (changing the rotation rate). The response is resonant when the denominator vanishes, i.e. when tj + mQ. — 0, which 
means that the tidal frequency in the inertial frame vanishes. 



A2 The case I = 3 

Uiaq 2 q s ojn 2 (R 2 



ai 



Wi 



3(ojiui 2 — 3q 2 Q 2 
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61 = 



14iaq 2 q 3 ujQ. 2 (R 2 - 2r 2 
3(wi£J2 — 3q 2 Q, 2 ) Rr 



C2 = 



l4aq 2 q 3 oj(ui + 2mD.)Q. 2 r(R 2 



3(wia;2 — 3q 2 Q 2 ) 



R 



W 3 = auifi 



2m 



56q 2 uiQ 



3(loilj 2 — 3q 2 Q, 2 ) 



14aq30JtJifi r 
9(cjia;2 -3q^ 2 )R' 

9R' 



The response is resonant when the common denominator vanishes, i.e. when 

2\ 1/2" 



-m ± 



9 -to' 



but not for |m| = 3, because then q 3 in the numerator also vanishes. 



(A3) 
(A4) 

(A5) 



A3 The case I = 4 

81iaq 3 q4UiujiQ 2 



02 



r(R 2 



9W1W2W3 - (27g|w 3 + 32q 2 uj 1 )Q 2 R? 

81aq 2 q 3 q40jn 3 (5R 2 - 7r 2 



62 



27iaq3g4CJOJif2 



(3i? 2 - 5r 2 



2i? 2 

„2 



W4 = aujQ, 



§u-ioj 2 uj s - (21q 2 m + 32g|wi)fi 2 2i? 2 

81aqiq 2 q 3 q4U}Q 4 r 2 (10_R 2 - 7r 2 ) 

~9ojiuj 2 uj 3 - (27g|cj 3 + 32q 2 u 1 )n 2 2R 2 

243gln(cjicj 2 - 3q%Q 2 ) 



C3 



9aiiW2W3 — (27g2^3 + 32g 2 a;i)fi 2 

81ag4a;r2(ix)iaj2 — 3g 2 ^ 2 ) 
$0JiUJ 2 m - (27q\u 3 + 32g 2 wi)n 2 &R 2 ' 

8laq 3 q4UjQ. 2 (ujl - 5<?ff2 2 ) r 2 (i? 2 



W 2 = 



9wia;2W3 - (27g§w 3 + 32q' 2 u) 1 )Q 2 



2i? 2 



2m - 



9wiW2a;3 - (27q\uj 3 + 32q 2 ui)n 2 



16R 2 



Again, the response is resonant when the common denominator vanishes, i.e. when ui satisfies the cubic equation 

42w 3 + 63mu; 2 Q, - 36(2 - m 2 )uj^l 2 - 4m(ll - 2m 2 )fi 3 = 0, 

but not for \m\ = 4, because then 94 in the numerator also vanishes. 



(A6) 
(A7) 
(A8) 
(A9) 

(A10) 



APPENDIX B: ANALYTICAL RESULTS FOR INTEGRATED RESPONSES 
Bl Homogeneous body with a solid core 

The equivalent of equation (|113[) for general values of I and m is 
r °° duj _ 3(21 + l)ne 2 [(l + l) 3 (l - l)(2l + 3)(l 2 - m 2 ) + l 3 (l + 2) (21 - l)((l + l) 2 



im[*r»]- 



l 3 (l + l) 3 (l - l) 2 (2l + 3) (21 - 1)(1 - q 2 ^ 1 ) 



(Bl) 
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Tides in rotating barotropic fluid bodies 



B2 Homogeneous fluid body with a solid core of a different density 

The equivalent result in the case I = m — 2 when the solid core has a density / _1 times that of the fluid envelope is 



T rf2, M du> 1007T 2 

Im# a (w) — = -— e 

w bo \ 1 — cr 



1+ 2 



We recover equation (|113[) either as / —¥ 1 or in the limit a 3 <C /. 



B3 Piecewise-homogeneous fluid body 

The equivalent result for the piecewise-homogeneous fluid model is 



dw 



IOOtt 2 



i^) (1 - /)2(1 - Q)4 ( 



1 + 2a + 3a + -a 



3 3^ 2 



1+1 ^-U 3 



For small / and small a, specifically a C / < 1, we again recover equation (|113 
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